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TWO  APPROXIMATIONS  OF  SOLUTIONS  OF 
KAMILTON-JACOBI  EQUATIONS 


H.  G.  Crandall  and  P.  L.  Lions 

Introduction.  The  main  results  of  this  paper  concern  the  approximation  of  solutions  of  the 
Cauchy  problem  for  first-order  partial  differential  aquations  of  Hamilton- Jacobi  type. 

Most  of  the  presentation  here  will  be  In  the  context  of  problems  of  the  form 


+  H(Du)  “0  in  ***  *  (0,») 


(IVP) 


u(x,0) 


Ug(x)  In  R 


where  H  e  C(R*^)  (the  continuous  functions  on  Uq  e  BUC(I^)  (the  bounded  and 

uniformly  continuous  functions  on  1^),  and  Du  -  ^'ix^ '  *  *  * spatial  qradlent 
of  u.  The  problem  (IVP)  Is  technically  simpler  than  the  "general  case"  In  which  the 
Hamiltonian  H  may  depend  on  x,  t  and  u  as  well  as  Du,  and  we  prefer  to  keep  the 
Ideas  clear  and  constants  simple  by  dealing  primarily  with  (IVP).  (See  the  comoents  In 
Section  4  regarding  more  general  equations.)  Two  sorts  of  approximations  of  (IVP)  will  be 
considered  here  -  finite  difference  schemes  and  the  method  of  vanishing  viscosity.  Before 
describing  these  approximations,  we  briefly  review  some  basic  facts  concerning  (IVP). 

Analysis  by  the  method  of  characteristics  shows  that  if  H  and  Uq  are  smooth  and 
Uq  Is  compactly  supported,  then  (IVP)  will  typically  have  a  unique  C  solution  u  on 
some  maximal  time  interval  0  <  t  <  T  for  which  lim  u(x,t)  exists  uniformly,  but  this 


t« 

limiting  function  Is  not  continuously  differentiable.  Thus  Du  "becomes  discontinuous" 
at  t  *■  T  (or  "shocks  form"),  if  one  Insists  upon  a  solution  of  (IVP)  which  Is  defined 
for  all  t  >  0,  It  Is  therefore  necessary  to  deal  with  functions  which  are  not  smooth.  On 
the  other  hand.  It  is  relatively  easy.  In  the  above  circumstances,  to  produce  Lipschlts 
continuous  functions  u  on  II  x  [0,w)  which  satisfy  (IVP)  If  the  equation  Is  understood 
In  the  "almost  everywhere"  sense.  However,  "generalised"  solutions  in  this  sense  are  not 
unique. 
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Recently,  a  way  of  identifying  a  uniquely  existing  solution  for  a  class  of  problems 
which  Include  (IVP)  as  a  special  case  was  given  by  the  authors  in  [2],  [3]  (see  also 
[1]).  The  relevant  solutions  of  scalar  nonlinear  first  order  equations  are  called 
"viscosity  solutions"  and  they  are  known  to  be  the  solutions  of  primary  interest  in  many 
areas  of  application  (e.g.,  optimisation,  control  theory,  differential  games  —  ).  See, 
e.g.,  [6],  (10],  [11].  The  term  "viscosity  solutions"  refers  to  the  fact  that  all 
solutions  obtained  via  the  method  of  vanishing  viscosity  are  in  this  class.  The  main 
properties  of  viscosity  solutions  relevant  for  the  current  work  are  recalled,  in  the 
context  of  (IVP),  in  Section  1. 

In  this  paper  we  will  approximate  the  viscosity  solution  of  (IVP)  l>y  solutions  of  the 
Vsseral  class  of  finite  difference  schemes  Introduced  below.  Indeed,  explicit  error 
estimates  are  given  relating  the  viscosity  solution  of  (IVP)  and  the  solutions  of  these 
finite  difference  approximations.  Ha  also  show,  under  suitable  hypotheses,  that  if 
c  >  0|  tt*  is  the  solution  of  the  problem 

f  1^  +  H(0u®)  -  eAu®  -  0  in  *"  X  (0,«)  , 

I  e  N 

V.  u®(x,0)  -  Uq(x)  in  r"  ; 

and  tt  is  the  viscosity  solution  of  (IVP),  then  |u*(x,t)  •  u(x,t)|  <  c  /~e  for  x  e 
and  t  >  0.  This  is  done  in  Section  5.  Estimates  like  this  have  also  been  obtained  in 
H.  H.  PlaaUng  [7]  and  P.  L.  Lions  [10]  by  indirect  arguments  involving  stochastic 
differential  games. 

He  now  describe  the  class  of  difference  approximations  to  be  considered  here.  For 
notatlonal  simplicity  only,  wo  will  assume  that  H  -  2  in  most  of  the  presentation.  The 
corresponding  definitions  and  results  for  general  N  will  be  clear  from  this  special  case, 
and  we  will  not  explicitly  formulate  them,  a  generic  point  in  R^  will  be  denoted  t)y 
(x,y)  and  we  will  write  Du  -  (u^.u^).  Given  mesh  sixes  Ax, Ay, At  >  0,  the  value  of  our 
numerical  approximation  at  (*j»yit.t„)  -  ( JAx,kAy,nAt)  (J,k,n  e  E)  will  be  denoted  by 
Capital  letters  U,  V,  ••••  will  denote  functions  on  the  x,y  lattice 
A  -  •  1'^  ®  their  values  at  (Sj.Yjj)  viH  be  written  • 
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Thus  o'*  r«pr«a*nts  tha  atata  of  our  numarical  approximation  at  tha  tima  laval  n^  and 
it  ia  a  function  on  A  with  valuaa  0?  .  .  Tha  notationa  X*  ■  At/Ax,  X^  «  At/Ay, 

J 


‘  "j.k  **• 


Tha  diacrata  approximationa  of  (IVP)  of  intaraat  hara  ara  axplicit  marchiny  achaaaa  of 
tha  form 

<11  -  ®‘""-p.)c-r . ' 

whara  p,  q,  t,  a  ara  fixad  nonnaqativa  intayara  and  G  ia  a  function  of  (p4'qt’2)(r*'a^2) 
variahlas«  (At  thia  ata9a  wa  ara  ignoring  tha  dapandanca  of  G  on  te,  ^  and  At.)  To 
aimplify  notation,  (1)  will  alao  ba  wrlttan  aa 

(2)  o"*’  -  i(o“)  . 

Wa  aay  that  (1)  haa  "dlffaranoad  form"  if  thara  axlata  a  function  g  auch  that 


®^**J-p,k-r'*  *  *  ,)t+a«-1 

(3) 

11  f]  fi  ^  fi 

_1  j+gfh+a+l  _♦  "j-p,)c-r  _♦  "j'*q<-1»W'a\ 

“’'■Ax  '•••'  Ax  »  Ay  tv  >  ' 

In  ordar  that  tha  achama  (1),  (3)  ba  conaiatant  with  tha  aquation  u^  ■•’  lI(u^,Uy)  >  0 
occurring  in  (XVP),  wa  muat  hawa 

(4)  g(a,...,ai  b,...,b)  •  H(a,b)  for  a,b  e  ■  . 

Whan  (3)  holda  wa  call  g  tha  nuamrical  Hamiltonian  of  tha  achaaw.  Finally,  wa  will  say 

that  (2)  (or  (3))  ia  monotona  on  (-R,R]  if  )  is  a 

nondacraaaing  function  of  aach  argument  aa  long  aa  Id'S),  I,  ,1  <  R  for 

J-p  4  14  J'fq,  )c-r  4  m  4  k'^atl,  J-p  4  !•  4  k-r  4  m*  4  k'*-s.  Roughly  spaaking,  R 

will  ba  a  priori  bound  on  |u^l>  Iwyl  for  the  solution  of  (IVP). 

Our  main  result  is 
2 

Thaoram  1i  Lat  H  i  >  ■»  R  ba  continuous  and  Ug  ba  bounded  and  Llpschits  continuous 
on  1?  with  L  as  a  Lipschits  constant.  For  x",  X^  >  0  and  fixad,  lat  the  schema  (2) 
have  diffaranead  form,  ba  monotona  on  (-(b-M  and  ba  consistent  with  (IVP).  Assuma 

tha  nuswrloal  Hamiltonian  g  is  locally  Llpschits  continuous.  Oaflna  u’  by 
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than  u",  n  •  1,2,...  by  (2).  I«t  u  ba  tha  vlacoaity  solution  of 

j#K  0  3  K 

(IVP).  Than  thara  Is  a  constant  c  depandlng  only  on  supIuqI,  L,  g  and  Ndt  such  that 

(5)  lu”  .  -  u(x  ,y  ,nfit)|  <  ci/Sb) 

j,K  3  a 

for  0  <  n  <  N  and  all  j,k. 


The  body  of  the  paper  Is  structured  In  the  following  wayt  Section  1  in  devoted  to  a 
review  of  tha  properties  of  viscosity  solutions  as  needed  herein.  Kxaaplea  of  difference 
scheaes  satisfying  the  assui^tiona  of  Theorem  1  are  presented  in  Section  2.  Theorem  1  is 
proved  in  Sections  3  and  4,  with  Section  3  consisting  of  preparatory  lewsis  on  the  mapping 
U  *  6(U)  while  Section  4  contains  tha  proof  of  <5).  Section  4  concludes  with  remarks  on 
variations  of  Theorem  1.  Tha  approxlMtlon  of  (IVP)  by  (IVP)  ^  is  treated  in  Section  S. 

We  bring  this  long  Introduction  to  a  close  with  some  resuirks:  First  of  all,  some 
convergence  results  are  given  in  s.  N.  Kruzkov  [8]  for  convex  Hamiltonians  H,  using  some 
estimates  available  only  in  this  special  case.  Next  (see,  e.g.,  P.  L.  Lions  [10]  Chapter 
16)  problasM  like  (IVP)  are  closely  related  to  hyperbolic  systems  satisfied  by  Du.  If 
N  •  1  this  relation  is  quite  simple,  since  if  u  solves  (IVP),  then  v  •  u^ 
solves  a  scalar  conservation  law: 

I  It  *  "ST  ° 

1  ’*'*0 

v(x,0)  -  v-(x)  »  (x)  in  K  . 

0  ox 

In  this  case  the  schemes  presented  here  are  related  to  those  studied  by  N.  G.  Crandall  and 
A.  Hajda  [4],  N.  N.  Kuznetsov  [8]  via  tha  corresponding  substitutloni 


v"  -  — 

'^k  Ox 


Vk 


u"  -u" 

k»1  k 
Ax 


Wo  refer  the  reader  interested  in  other  aspects  of  HamiltonrJacobi  equations  to  M.  G. 
Crandall  and  P.  L.  Lions  [2],  P.  L.  Lions  [10],  and  M.  G.  Crandall,  L.  C.  Evans  and  P.  L. 
Lions  (1).  Finally  we  reawrk  that,  in  an  ongoing  investigation,  P.  Souganidls  [11]  has 
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WP>w>wm,., 


formulated  ganaral  aroroxlmatlon  raaults  which  appaar  to  apply  alika  to  dlmanalonal 
aplittlng,  Mic-mln  representatlona,  approximation  by  (IVP)^  and  numerical  achamaa.  Hia 
argumanta  ara  ralatad  to  thoao  given  harain* 
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S«ctlon  1.  Viscosity  solutions  of  (IVP). 

As  rsesllsd  in  ths  Introduction,  one  cannot  solve  (IVP)  in  a  classical  way  on 
***  X  (0,»)  in  general,  while  Lipschitr  continuous  “generalised  solutions*  in  the  alnost 
everywhere  sense  exist  under  mild  aasusiptions  but  are  not  unique  (exanples  are  given,  for 
instance,  in  C2]).  The  resolution  of  these  difficulties  is  given  in  [2],  the  results  of 
which  ioply,  in  a  round-about  way,  the  theorem  stated  below.  This  theorem  is  proved 
directly  in  [1].  (We  use  x  to  denote  points  in  below.) 


Theorem  (Ixistence  and  uniqueness).  Let  H  6  C(«*‘),  Uq  e  BUC(I^).  Then  there  is  exactly 
one  function  u  9  »  (0,T] )  for  all  T  >  0  such  that  u(x,0)  -  Uq(x)  and  for  every 

♦  e  c’(B^  X  (0,»))  and  T  >  Oi 


(1.1) 


If  (Xg,t||)  la  a  local  maximum  point  of  u  -  # 
\  on  X  (0,T],  then 


and 


(1.2) 


If  A  local  sdnimum  point  of 

<  u  -  ♦  on  l/*  X  (0,T),  then 

U  (x„,tg)  +  H(Df(Xg,t„))  >  0  . 


The  function  u  whose  existence  and  uniqueness  is  asserted  by  the  theorem  la  called 
the  viscosity  solution  of  (IVP).  A  continuous  function  u  on  x  [0,T]  which  satisfies 
(1.1),  (1.2)  is  called  a  viscosity  solution  of  t)ie  equation  u^  H(Du)  -  0  on 
1^  X  [0,T].  See  [1],  [21  for  the  appropriate  notions  for  more  general  equations.  There 
are  also  useful  equivalent  ways  to  fomnilate  the  notion  of  viscosity  solutions  ( [  1  ] , 

[2]).  JMong  the  desirable  properties  of  the  notion  of  viscosity  solutions  is  its 
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consistency  with  the  clssslcsl  concept.  Thst  Is.  if  u  is  a  classical  (i.e.,  C  ) 
solution  of  u^  H(Du)  •  0,  than  it  is  a  viscosity  solution  and  if  u  is  a  viscosity 
solution  than  u^(xg.tg)  *  H(Du(xg,tg))  •  0  at  any  point  (Xg.tg)  where  u  is 
differentiable. 

The  other  Infonaatlon  we  want  to  recall  consists  of  various  estimates  on  the  behavior 

of  solutions  of  (IVP).  To  record  these,  for  each  t  >  0  let  S(t)  i  BUC(B^)  *  be 

the  time  t  map  associated  with  (IVP).  That  is,  S(t)uQ(x)  •  u(x,t)  where  u  is  the 

viscosity  solution  of  (IVP).  ffe  also  put  ifl  ~  sup|f(x)|  and  f'*'  -  max(f.O).  The  next 

E** 

result  follows  from  [21i  see  also  [1]. 

Proposition  1.1.  Let  H  e  C(E^)  and  S(t)  be  as  abovei  Ug.Vg  e  BUC(E^).  and  t  >  0. 
Than 

(i)  l(S(t)Ujj  -  S(t)V^)‘^l  <  l(Up-Vjj)*l  . 

(il)  IS(t)U||  -  S(t)V||l  <  lUg-Vgl  . 

(ill)  inf  Ug  <  tH(0)  S(t)u.  <  sup  u.  . 

if  e“ 

(iv)  |8(t)Ug(x4y)  -  8(t)Ug(x)|  <  sup  |Ug(t+y)  -  Ug(s)|,  for  y  e  e"  . 

seE" 

(v>  If  Ir  is  a  Ilipschits  constant  for  Ug.  then  it  is  also  a  Upschltc  constant  for 
8(t)Ug  and  l8(t)u.  -  8(T)Ugl  <  (t-T|  sup(tH(p)|  i  |p|  <  b)  . 

e" 

Tht.  key  point  here  is  (i).  The  astisMte  (1)  iRiplles  (11)  upon  using  (i)  with  ug 
and  Vg  Interchanged.  Clearly  (i)  also  implies  S(t)U|j  >  S(t)v^  if  u^j  >  v^,  which  In 
turn  implies  (iii)  since  v  •  c  -  til(O)  is  a  classical  (and  so  the  viscosity) solution  of 
(IVP)  with  the  constant  initial  datum  c.  Choosing  Vg  •  sup  Ug  or  inf  Ug  and  using  the 
order  preserving  property  yields  (ill).  Next  (iv)  follows  from  (li)  Isecause  S(t)UQ(  M-y) 
is  the  solution  of  (IVP)  for  the  initial  datum  Ug(*'^).  Since  (iv)  shows  that  a  modulus 
of  continuity  for  Ug  is  also  a  modulus  for  S(t)Ug(«),  the  flri:  .iicrtlon  of  (v)  is 
clear.  The  Lipschits  property  in  the  time  is  easily  deduced  from  ths  !'(>.<  t  ion 
*  H(Du)  '0  in  the  viscosity  sense  (see  (2))  or  in  other  ways. 


S«ctlon  2»  BxaiiVlea 


Ha  ba^in  with  N 

(2.1) 


1  and  wrlta  (IVP)  In  tha  form 


f'  M  *  H(u  )  “  0,  for  t  >  0,  X  e  R  . 

X 

^  u(x,o)  “  X  e  m  , 


in  thia  casa.  An  the  first  axampla,  wa  considar  the  schema 


whara  6  >  0  is  givan.  Tha  ralation  (2.2)  piay  be  rewritten  as 


Bailing  the  differenced  form  clear.  The  numerical  Hamiltonian  is  given  by 
g(a,0)  «  H((at«'S)/2)  -  (0-o)6/A*  for  a,fi  e  R.  clearly  g(a,a)  "  B(a)  and  so  (2.2)  is 
consistent.  The  acheaa(2.2)  is  monotone  on  [-R,R]  if  1  -  26  >  0  (monotonicity  in 
0^),  and  8  -  l*|H'(a)l/2  >  0  for  |a|  <  R  (monotonicity  in 

relations  are  achieved  by  first  choosing  0  <  6  <  1/2  and  then  sufficiently  small. 
This  schesw  is  analogous  to  the  XiSx-Prledrlchs  scheme  for  conservation  laws,  see  [41. 

In  a  similar  way,  the  schemas 

(2.3)  -  AtH 
or 

(2.4)  0^’  -  0^  -  AtB 


liave  the  desired  properties  if  H  is  nonincraasing  (for  (2.3))  if  H  is  nondecreasing 
(for  (2.4))  and  1  >  l’‘|H'(a)|  for  |a|  <  R.  Those  are  simple  "upwind*  schemas.  Next, 
let  e  R  and  assume  H'(aKa>a|j)  >0,  so  R'  changes  sign  at  a  •  c^.  Set  8(s)  > 
if  s  <  Cig,  6(s)  •  0  if  s  >  Ug. 


He  consider  the  schi 


r 


f 


t. 


The  numerical  Hamiltonian  la  now 

g(o,e)  -  e(B){H(B)  -  H(ag)}  +  (1-0( o) ) (h( a)  -  H(o^))  ♦  0(0^)  » 

and  thus  (2.5)  is  consistant.  Remarking  that  g(a, B)  may  be  written  as 

g(a,B)  •  H(B  *  0^)  H(a  v  a^)  -  HCag),  where  "a"  and  "v*  denote  "maximum*  and 

"minimum",  It  Is  clear  the  g  Is  locally  Llpschlts  if  H  is  locally  Llpscbitz.  ylnally 

one  checks  that  (2.5)  la  sunotone  on  (-R,'tR]  If  1  -  l’‘|H'(a)|  >  0  for  |  a|  <  R. 

In  fact,  all  the  above  examples  are  merely  adaptations  to  Hamllton-Jacobl  equations  of 

well-known  schemes  for  conservation  laws  via  the  remarks  in  the  Introduction.  As  explained 

in,  e.g.,  M.  G.  Crandall  and  A.  Majda  [4],  there  is  a  class  of  schemes  for  the  conseirvation 

law  +  (H(v))jj  “  0  called  "monotone,  in  conservation  form"  with  the  following 

structure;  ”  v"  -  X*4’*g(v”  )  -  now  the  function  g  is  called  the 

3  3+  3-p  3+q 

"numerical  flux*.  Consistency  then  means,  as  before,  g(a,...,a)  -  H(a),  and  monotonicity 
means  that  the  map  (v"  ►*  ^  is  a  non-decreasing  function  of  each  V^,  Then  we  may 

write  the  corresponding  scheme  for  the  approximation  of  (IVP), 


”i  . ■ 


which  is  in  differenced  form  and  consistent. 

Next  if  N  >  1,  the  relation  between  (ivp)  and  conservation  laws  disappears.  For 
N  «  2  (to  simplify,  as  always)  we  mention,  the  analogue  of  the  Lax-Frledrlchs  scheme,  i.e. 


(2.6) 


f.hw 


V 


u"  u"  -  II" 

3-*-1.k  ”j-1.k  “j,k-»-1  ^1.k-1 

2Ax  '  2iy 


-  g-  (  Vt,k  -  -j,k)  e  I  “j,k>i  -  “j.k-1 


-  20 


20 


Ax 


Ay 
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which  is  consistent,  in  differenced  form.  It  is  also  monotone  on  [-R,R]  provided  that: 

0  <  e  1  -  X*|h'(o,8)|  >  0,  and  0  -  X^|H‘(a,8)|  >  O  for  |  a| ,  |  B|  <R,  where  H'. 

4  1  2  ^ 

denotes  the  derivative  of  H  in  its  i^**  argument. 

As  a  last  example,  let 


i.k 


.  v" 

'  j+q+1 ,k+s+1 ■ 


j.k 


aX 

-  At  g,  (-LJbfiaizE 


Ax 


A^'v" 

+  1+a.k+s+1 

.  ^ - » 


aV 

Ay 


*  i+q-t-l  ,kts  \ 
Ay  ' 


and 


w"*’  -  G  (w"  w"  » 

J,k  2*  j-p,k-r'***''*j+<j+l.k+a+1 


A*W^  A^w” 

. 


*  i-t-q-n  ,k-*-s  •> 
'  Ay  J 


be  differenced  form,  monotone  on  [-R,R]  schemes  consistent  with  v  +  -  H  (v  ,v  )  -  0 
,  t  a  1  X  y 

“t  7^  ”2*’'x'''y’  “  °  r«8P«ctively,  where  oe  (0,1).  Then  the  scheme 


0^’  -  a  G  (u"  . 

J»k  1  j-p,k-r 


j+<It1  ,k+s+1 


)  +  (1-o)G,(o' 


2"^J>p,k-r' 


'°j+q+1,k*e+1* 


differenced  form,  is  monotone  on  (-R,R]  and  is  consistent  withi  +  H  (u  .u  )  + 

Ot  1  X  y 

”2^“x'“v*  ”  we  could  build  schemes  in  this  way  for  +  H,  (-Ir)  t  H  C-^)  =  0 

^  1  OX  2  9y 

from  schemes  for  |S  *  (|H,  .  „  ^  ^ 
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Section  3.  Stability  Properties  of  the  Sehaaea 

The  notation  In  this  Section  aasuJ&es  two  space  dlnenslons>  but  everything  herein 
easily  generalises  to  arbitrary  dimensions.  By  capital  letters  U,  V,  etc.,  we  denote 
bounded  bllnflnlte  double  sequences  {U.  }  ,,  l.e.  bounded  functions  on  the 

m  2 

set  of  all  such  being  i  (S  )  which  we  equip  with  the  norm 

lul  -  sup  |U.  .  I  . 

Let  G  be  the  self-Mp  of  t  (S^)  defined  by  (3),  l.e. 


(3.1)  G(U) 


J»h  J»k 


u  -  ^"j-P>lt-r . ^Vq>1,k-t-s^ 


Wa  now  invaatlgata  tha  propartiaa  of  G*  By  daflnition  (2),  (3)  la  monotone  on 

(-R,Rl  if  the  restriction  of  G  to  C  -  {o  e  «  |A*U,  |  <  Rte,  1  <  R/^, 

for  t,a  e  s}  preserves  the  natural  ordering  of  i  (S  ).  We  Identify  1  e  R  with  the 

constant  function  X  on  S^.  From  the  fora  of  6,  It  Is  clear  that 

G(0  +  A)  -  5(0)  *  A  for  all  o  8  1*(»*)  and  A  e  R  , 

that  Is  G  cosaiutes  with  tha  addition  of  constants.  Now,  It  Is  a  siaple  fact  that  order- 

2 

preserving  aapplngs  coaautlng  with  the  addition  of  constants  are  non-expanslve  In  t  (S  ) 
(N.  G.  Crandall  and  L.  Tartar  [5]).  Indeed,  If  0,  V  e  C  then  u  <  V  -f  A  with 
A  >  l(u-V)^l_.  But  V  4'  A  e  C  and  thus  (using  aonotonlclty)  we  deduces 

G(U)  <  5(V  +  A)  -  G(V)  +  A,  so 

(3.2) 

I{G(0)  -  G(V)>  1^  <  l(0-V)*l_ 


which  lapllas  In  particular,  the  non-expanslveness  on  C. 

Another  siaple  property  of  G  la  that  it  cooautes  with  translations,  l.e.  if 


l,a  e  X  are  fixed  and  t  Is  the  linear  mapping  defined  by:  ( T.  U). 

*,a  t,m  j,K 


0, 


then  we  have  t  g(U) 
i,a 

with  (3.2)  we  can  deduce 


G(t 


t,B 


J+t,k+m 

U)  for  all  vet  (X^).  Using  this  property  in  conjunction 
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t. 


(3«3) 


(IA“  G(U)I,  <  lA*  Ul-  , 
lA^  G(0)  1^  <  lij  01^  , 


for  all  U  e  C<  For  example, 

lA*  G(U)I_  -  It,  „G(U)  -  G(a)l  -  IG(T  U)  -  0(0)1 


<  IT^  gO-Ol,  -  ia“oI„  . 

since  T,  leaves  Invariant  C  for  any  t,m  e  S.  An  Inaaedlate  consequence  of  (3.3)  is 
t,m 

that  G  leaves  C  invariant  (G(c)  ^  C). 

The  last  property  we  observe  is  the  followinqi  If  n,j  >  0  and  0  6  C,  then 


Ig"*^(0)  -  g''(0)I^  <  IG^(O)  -  01^ 


<  IG^"‘'(0)  -  G^”‘**’’(0)  <  j  IG(U)  -  Dl„ 

t-0 


and  from  the  explicit  form  (3),  we  finally  obtain; 

(3.4)  Ig”*’^(0)  -  g"(0)I„  <  j  At  K  for  n,j  >  0,  0  e  C 
where  K  is  given  by; 

(3.5)  K  -  sup{|g(e,  „)Ii  0  <  t  <  ptq+1,  0  «  m  <  r+s+1,  15,  _|  <  R)  . 
We  record  all  these  properties  in 


Proposition  3.1;  Let  R  >  0  and  c  be  defined  as  above.  Let  the  scheme  be  sonotone  on 
♦  *  2 

[-R,R]  and  G  ;  C  >  t  (Z  )  be  given  by  (3)  with  g  bounded  on  bounded  sets.  Then  we 


have; 


(1) 

G(U)  < 

G(V)  for 

u,v  e 

C,  0  <  V  ; 

(ii) 

G(U+X) 

»  G(U)  + 

X  for 

0  e  c,  X  e  z  1 

(iii) 

■G(0)  - 

-  G(V)I„  < 

IU-Vl_ 

for  u,v  e  C  1 

(iv) 

lA*  G(0)l,  <  lA* 

ul,,  lA^  G(0)l„  <  IA][ 

(V) 

G(C)  e 

C  » 
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(Vi)  Ig'**^(U)  -  g"(U)I,  <  j  At  K  for  0  6  C,  n.J  >  0 
where  K  Is  given  by  (3.5) i 
(vii)  Ig"(U)I^  <  lUl,  +  n  At  A,  tor  U  6  C 
where  A-  |g(0,...,0)|  -  |H(0,0)|. 

The  only  new  property  is  (vii)  which  is  essily  proved  observing  first  thst 
q"(0)  •  -n  At  g(0,0,...0)  end  so 

IG®(0)I,  <  Ig“{0)  -  g"(0)I,  +  Ig"(0)I, 

<  lOl.  +  n  At  A  . 
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section  4.  Proof  of  theorea  1 

The  proof  of  Ttieorea  1  given  here  is  related  to  the  proof  of  uniqueness  of  viscosity 
solutions  of  (IVP)  presented  in  [1]and  it  also  involves  estimates  introduced  in  [2]. 
Throughout  this  section  we  will  assume  the  hypotheses  and  notation  of  nteorem  1 .  In 
addition,  we  will  at  first  assume  that 


(4.1) 


{ 


u(x,y,t) 


as 

as 


|x|  +  lyl  ♦  •  and 
Ijl  +  |k|  *  • 


hold  uniformly  for  bounded  t,  nAt  >  0.  This  assusqption  allows  a  simplified  presentation 
and  la  easily  relaxed  later.  Moreover,  (4.1)  holds  if  u^  0  as  |x|  +  |y|  *  *. 

It  will  t>e  convenient  to  define 
(4.2)  a  -  *  (0,-),  Qj  -  X  tO,Tj 

and  the  discrete  analogues 


(4.3) 


“  A  X  (At«*)  «  {(x^,yj^,nAt)  :  j,)c  -  0,  ±  1,...>  n  •  0,1,...,} 
“  {(Xj,yj^,nAt)  e  a**  j  n  <  n)  . 


Hereafter  T  >  0  and  N  e  Z  satisfying 

(4.4)  (N-1 )At  <  T  <  NAt 

are  fixed.  We  seelt  to  estiute  u(x,,y  ,nAt)  ~  U?  ^  this  end,  we  will  assume 

J  X  j*x 


sup  (u(x  ,y  ,nAt)  -  O.  .  )  •  O  >  0 


(4.5) 

0<n<N 

and  then  produce  an  upper  bound  on  c.  In  exactly  the  same  way,  if 

inf(u(x  ,y  ,nAt)  -  D?  )  "  -o  <  0,  we  could  estimate  o  and  the  conjunction  of  these 
3  K  ],K 

estimates  Isounda  (u(x  ,y  ,nAt)  ~  U4  .  )• 

J  X  J,K 

We  are  going  to  define  a  function  <1*  <  B  x  a  which  is  a  principal  ingredient  in 

the  proof.  This  function  depends  on  a  and  T  above  as  well  as 

(4.6)  M-  sup  (lu  (5))  tTlHlO.O)!)  *  1  , 

ieie«* 


a  positive  parameter  e,  and  a  function  S 


R  X  R  «’  R  as  follows:  For 


(  C,t)  e  a  and 
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(n,s)  e  g 


-  u(C,t)  -  a"  ^  (t+«)  +  (5M  ♦  |)6g(C-n,t-«) 
whar*  (nfS)  -  (x^.Yj^.nAt)  and  -  8(C/e.t/e)  . 


Wa  raaurk  that  (4»6)  guarantaaa 


(4.8)  |tt|  <  M  on  g^  and  <  H  for  0  <  n  <  M  i 

aaa  Proposition  1.1  (iii)  and  Propoaition  3.1  (vii).  Tha  function  B  of  (4.7)  will 

satiafy 


B  la  aaooth  on  R  x  r,  o  <  8  <  1,  8(0,0)  •  1 

(4.9) 

and  8(€,t)  -  0  it  U|^  ♦  t^  >  1  , 

aa  wall  aa  othar  oonditiona  iapoaad  latar.  tha  naxt  atap  la  to  ■axiaiaa  4*  ovar 
®T,8  “  ®T  *  ®M* 

Lawa  4.1.  Ondar  tha  abova  aaauaptiona  thara  ia  a  point  j,  auch  that 

<i)  >  t(C«t,n.a)  for  (C.t.n.a)  e 


(11)  B,(Co  -  hj.  to  -  aj,)  >  3/5  . 

Proof.  Tha  axlatanca  of  ((Q>to,no>ao)  followa  froa  tha  fact  that  If 
(5*/t*,n*,a*)  «  (Ij  „  «nd 

(4.10)  t(5*,t*,n*,a*)  ♦  aup  ^  , 


than  (C^/t*,i)^,a^)  raaaina  boundad.  Indaad,  froai  (4.5),  (4.7),  (4.9)  ona  aaaa  that 


(4.11) 


aup  ♦  >  aup  (u(x.,y  ,ndt)  -  a"  )  -f  5M  -  a  ♦  5M  , 
$i-„  J.kCR  ^ 

0<n<H 


whila  (4.7),  (4.8)  laply 


T(e,t,n,a)  <  211  if  B,(t-n,t-a)  -  0  . 


Hanca  (4.10)  lapliaa  Bo(C*“n*,t*-a*)  >  0  for  larga  t  and  ao  (?*-n*)*  (t*-8*)^  <  e* 

Proa  thla,  (4.1)  and  (4.7)  wa  aaa  that  if  *  |n*|  *  *,  then 

lia  aup  T(C*»t*,n*,a*)  <  5H  *  c/2.  Thla  contradicta  (4.10),  (4.11). 


To  prova  (11),  observe  that 

2M  +  (5M  +  <»/2)6e(eo-'Vl'V*o’  ^  W5o'’^0''\)'‘o’ 

>  sup  f  >  5m  +  o  , 

so  Be<to"'’o'^o"“o’  *  *®'*  *  2o)/(10M  +  o)  >  3/5. 

In  what  follows  ve  will  put 

(4.12)  e  -  (Ax  +  Ay  +  At)’^^  -  (A*  +  +  1) ’'^^( At) , 

although  we  will  not  use  this  relation  except  at  certain  points  In  the  arguaent. 

There  are  now  several  cases  to  be  considered.  These  are:  tQ,Sg  >  0i  tg  >  0,  Sg  ~ 

0:  and  tg  0,  Sg  >  0.  We  begin  with  the  case  tg,Sg  >  0. 

1st  case;  tg  >  0,  Sg  >  0. 

It  follows  frost  lasBut  4.1  (1)  that  (Cg,tg)  Is  a  maxinwa  point  on  of 

(C,t)  ♦  u(C,t)  -  ^  t  +  (5M  +  ■j)8g(t-ng,t-aj|).  By  the  definition  of  viscosity  solutions  and 
the  aquation  solved  by  u,  we  therefore  have 

f  -  (5M  +  1)0^36* * 

(4.13) 

H(-(5M  +f)D58^(5g-»\,.VSp))  <  0  . 

(Mere  and  tielow  D,.  6  «  0,6  stand  for  the  Indicated  derivatives  of  P  (C,t)  which  are 
then  evaluated  at  the  point  shown.) 


The  analogous  estlnate  on  the  discrete  side  requires  eore  work.  Let 


(4.14) 

(hg.Sg)  -  ((X  ),  ngAt)  i 

•*0  0 

then 

minimizes  (j,k,n)  ♦  +  ■|(nAt)  -  (5M  + 

over 

nAt) 

e  Q„.  Thus 

N 

(4.15) 

u" 

J.k 

>  “v»,  •  J  ‘v«'“ 

-  (5M  ♦  f><Bj(Cg-ng,t^-Sj,)  -  8^(eg  -  (JAx.kAy),  tg  -  nAt)) 
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for  0  <  n  <  N  and  all  j,k.  We  next  uae  the  Bonotonlclty  of  G  and  (4.15)  with 
n  “  njj  -1  to  conclude  that 


"o  ♦  "o"’ 

"l  k  “  N  k 

Jo'*o  ^O'^O 


(4.16) 


*  "Vko  *  !  -  (5M  t  f)(B,<eo-n„.V»0>  -  V-0  " 

-  Bt  ,(<5Hy-)  .  ((j^.p,Ax.  (kg-DAy),  -  -0  >  At)  . 


•••»  V'o  *  *<^>'***>  • 

To  guarantee  the  validity  of  this  step  we  must  show  that  the  arguments  of  g  above  lie 
in  [- (L>1 ) ,Lt1 1 ,  since  the  monotonicity  is  only  assumed  in  this  case.  Consider  a  typical 
argument  of  g  above;  e.g. ,  any  of  the  x-differences  )uis  the  form 
(4.1?)5J*^(6^(?,,-nQ  -  ((ttl)Ax,mAy),  t^j-Sg  +  At)  -  Bg(?g->\j  -( lAx.mAy) .t^-Sj  +  At)) 
where  t,  m  are  bounded  Integers, 

(4.18)  1  +  |t|,  |m|  <  max(p,r,<i*’1,st1 )  +  1  -  IC  . 

Clearly,  the  difference  between  (4.17)  and  -(5ltto/2)D^ Bg( tQ“>\,,tp-Bu) ,  where  D,  denotes 

differentiation  in  the  first  spatial  argument,  is  estimated  in  t)ie  form 
-2 

cons,  e  (AxtAytAt),  where  the  constant  involves  bounds  on  the  second  derivative  of  0 
but  is  independent  of  e.  Invoking  (4.12)  we  find  these  errors  to  be  at  most 

(4.19)  cons.  (At)’^^ 

so  (4.17)  differs  from  -(5M  ♦  BglCg-hg,  tji-Sp)  by  at  most  1  if  At  is 

sufficiently  small.  Using  the  next  lemma  and  these  remarks,  we  see  that  the  arguments  of 
g  in  (4.16)  lie  in  [-(L>1),Lt1]  if  At  is  sufficiently  small.  (Parts  (li)  and  (Hi)  of 
the  lemoM  are  used  later.) 

tesBsa  4.2.  Let  (Co'^g''lg,ag)  be  as  in  Lemsu  4.1  and  L  be  the  Lipschlts  constant  of 
Ug.  Than 

(i)  (5lH-O/2)|Dj0^(eg-nj,t,,-Sg)|  <  L  . 
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Lat  Lf  •  Mx(|H(p)|  I  IpI  <  L)  and  tg  >  0.  Then 

(ii)  -  (5M>a/2)o^8^(Cg-ng.v»o»  <  ^1  -  ^  • 

If  also  T  >  tg,  then 

(ill)  (5H+0/2)|D^eg(5g-ng,tg-ag)|  ^  * 

Proof!  By  assunptlon,  the  mapping  5  *  u((,tg)  *  (5M+o/2> E-hg/tg-Sg )  is  maximized  at 

2 

C  *  Thus  for  C  e  R 

(5M+o/2)(6(E-ng,tg-Sg)  "  6 ( Eg ' hg , tg -8g ) )  <  u( Eg  , tg  )  -  u(E,tg) 

<  LUg-^l  . 

where  the  last  inequality  is  from  Proposition  1.1  (v).  The  inequality  (i)  follows  at 
once.  Similarly,  t  ♦  u(Eg,t)  “  ^  ^  <5M+o/2)$^(Eg~»^ ,t-8g )  is  maximized  over  [0,T] 

at  tg  >  0,  so  for  small  h  >  0 

(511+0/2)  [B^(Eg-ng,tg-h-Sg)  "  t  'o  ’  'Ip  ^  ^  ^  "O  ^  ^ 

<  (U(Eg.tg)  -  U(Eg.tg-h))  '  ^ 


Where  the  last  Inequality  is  from  proposition  1.1  (v).  The  inequality  (ii)  follows  at 
once.  If  also  T  >  tg,  then  one  malces  the  two-sided  estimate  in  the  obvious  way. 

Mow  tre  return  to  (4.16).  In  this  expression  we  replace  each  difference  in  the 
arguments  of  g  by  the  corresponding  derivative  of  -(5M+0/2)  B^( E-^^ ,  ^g~*g)  ^ 

thereby  creating  errors  we  can  estisMte  -  using  the  locally  Lipschitz  property  of  g  -  by 

•2 

multiples  of  (dx+dy+dt)e  .  Then  we  use  the  consistency  of  the  numerical  Hamiltonian  g 
with  B  to  conclude  that 


^  <  (511+0/2) 


At 


+  H(-(5M+O/2)D^0^(Eg-ng.tg-Sg)) 


^  c(AstAr.4t, 
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«h«r«  C  ia  *  conatant  ona  could  aaally  aatlaata.  Making  alMllar  arguMnta  on  tha  t- 
dlfCaranca  abova,  aa  furthar  daduca  that 

f  «  -(5MtO/2)D^B,(  V\)'S-»0> 

(4.20)  H(-(5Mta/2)D^e^(Cg-i\,,tj-aQ)) 

>  e(^y-S)  . 

Takan  togathar.  (4.13)  and  (4.20)  ylald 

„  <  C(tetyt) 

(with  a  naw  conatant  C).  Na  again  Invoka  (4.12)  ao  that  thla  bacoaaa 

a  <  C(Ax+Ay+At)''2-  C(X*  ♦  ♦  I)'''*  (At)'''* 

which  aatabllaliaa  tha  daalrad  aatlawta. 

Wa  turn  to  tha  caaaa  In  idilch  ona  of  tg  or  a^  la  O.  In  thaaa  eaaaa  wa  do  not 
naad  to  uaa  tha  Infonwtloa  that  u  la  a  aolutlon  or  tlta  datallad  propartlaa  of  6.  Ma 
raly  an  (4.11)  and  ali^ia  conaldaratlona  of  continuity.  Bowawar,  wa  will  raatrlct  B  to 
aatlafy 

(  B(C,t)  -  1  -  (lU*  ♦  t*)  tot  + 

(4.21)  \  and 

I  B(5,t)  <  i  for  Kl*  ♦  t*  >  -I  . 


In  tlia  avant  (4.21)  holda  wa  know  fro*  Laanut  4.1  (11)  tliat  *  (t^-ai^)^  <  ^  c^. 

Ranoa 


°5*e**o"'’o'*o"*0*  "  "  ^^o'^o’ 
Dt^i'^o'^O'^o'^O*  “  "  *S""o’ 


and  than  from  tmmm  4.2  wa  eoncluda  thati 

(4.22)  *  **  *-'*^®***'®’  » 

If  t„  -  T,  than  It.-a.l  <  a*(l,  -  ^)/(10llto)  , 


(4.23) 


(4.24) 


If  T  >  tg  >  0,  than  <  e^(L^  ♦  ^)/(10H+O)  . 

2nd  caaai  tg  >  0,  Sq  ~  0. 

By  (4.11).  Proposition  1.1  (v)  and  the  choice  of  0^ 

5M+0  <  ♦(Co.tjj.n^.o)  <  lultg.tjj)  -  u(i\j,tj)|  ♦  ludy^.tp)  -  u(iyj,o)| 

♦  (5lHO/2)B^(«g-.^.t„-.^) 

<  Lieg-npl  +  L^tji  +  (5lH-<)r/2)  . 

The  estlMte  (4.22)  holds  and  either  tg  -  0  or  one  of  (4.23).  (4.24)  holds.  The  al>ove 

2 

thus  yields  o  <  const,  e  .  and  again  wo  are  dona. 

3rd  case I  tg  ■  0.  Sg  >  0. 

in  a  Banner  slBllar  to  the  above  we  find  (using  also  Proposition  3.1  (v))  that 

(5IW-0)  <  ♦(tg.O.ng.Sg)  <  ^  |  {g  >  |  ^  ^Sg  ♦  (  5|«-  0/2  ) 

and  so 

O  <  2I-IV’\)'  *  ^>^0  • 

He  now  need  to  estlBate  Sg  suitably  and  Invoice  (4.22)  once  aore  to  coiq>lete  the  proof. 
Frocs 

♦  (Cg.tg.hg.Sg)  <*  ^  (  Cg  .  0  .  .  Sg  )  P  W  5g  .  0  .  .  Sg  “  At  ) 

we  conclude  that 

■  ‘5'^«/2)6^(«g-.y,,.Sg)  > 


where  the  notation  (4.14)  Is  being  used.  Since  8g( ^-lyj.tg-Sg)  -  '”"o*  * 

Lobsm  4.1  and  (4.12).  (4.21)  hold,  we  can  assuae  that  0^  has  the  quadratic  fom  (4.12) 
by  taking  dt  saall.  The  above  then  becosws 


,  2  .  ,  „2,  .  „V’  .  ^  -  i  „  <  . 

,200  In'^n  In'^n 
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Analysing  this  inequalityt  ws  concluda  that  Sg  <  contt«(c  *  tt).  Using  this  in  th« 

2  ^/i 

estimate  on  a  above  yields  a  <  const,  (c  *  At)  and  again  e  ■  (Ax+Ay^-dt)  *  yields 
o  <  const.  (Ax+Ay+At)'^  . 

Proof  of  the  General  Case 

The  remaining  step  la  to  remove  the  restriction  (4.1).  There  are  two  possible  ways  to 
do  this.  He  may.  for  example,  follow  the  unigueness  proof  in  [1)  and  replace  ^  by 

♦  -  ♦  +  2fi4(C,t,n,B) 

where  5  >  0,  C  -  ?(5+C<,/tttjj,y+y^,s+eg)  and  C  €  C*(a*>ill  x  R^xr)  ,  0  <  C  <  1i  C(0)  -  1 
and  *  point  such  that 

Then  adapting  the  above  proof,  one  reaches  the  desired  conclusion,  (see  also  (111). 

Another  argument  suikas  use  of  the  hyperbolic  nature  of  the  problem,  naMly  the  finite  speed 
of  propagation.  Observe  that,  without  loss  of  generality,  we  may  assume  R(0)  «  0 
(replace  u(x,t)  by  u(x,t)  >  tH(0)(  H  by  B  •  R(0),  g  by  g  -  H(0),  U*'  by  U*'  nAt). 

In  the  statement  of  the  next  result,  which  was  proved  in  [2] ,  we  use  the  notation  R(x,R) 
for  the  closed  ball  in  rN  with  center  x  and  radius  R  and  put  >  B(0,R). 

Theorem;  Let  ug,  VgCH'  (R),  let  H  6  (B  )  and  denote  the  semigroup  solving 

(ivp)  by  S(t).  Then,  if  yg  e  and  Ug(y>  •  Vg(y)  on  B(yg,R)  we  have 

(8(t)Ug)(x)  -  (S(t)Vg)(x)  on  B(yg,H-vt) 

where  v  -  IH' I  ^  and  r  ■  Hax(lDu.  ■  •Ov.l  Jl. 

L  (B^)  L  L 

On  the  other  hand,  it  is  clear  from  the  definitions  that  if  ~  k 

lj-jgl  <  R*  ll-jgl  <  B  (R  e  «)  then  0^  for  lJ-Jg|  <  R-Kn,  |k-k^|  <  R-Rn 

with  K  ■  max(p,r,<{ti,r*'1>.  Thus  if  Ug  S  v^^  in  B(Cg,R)  and  Cg  -  **0’^0*'  *"** 

v”  are  the  discrete  approximations  generated  by  our  scheme,  we  see  that  u|?  .  -  v" 

J#ic  Jpic 

for  jp  k  satisfyin^s 


f. 


*0  R  ^0  R 

II  -  s;' " 2  I’' - ‘ -5? 


2  -  Kn 


!••• 


IJAjc-x  I  <  (R-2Ax)  -  5-  (i»At), 
u 


IkAy-y^jl  <  (R-2Ay) 


-  —  (nAt). 


It  is  then  easy  to  conclude  the  arguiaent  by  roaarlclng  that,  uniformly  in  Zg  e  R^,  we 
may  find  R  large  enou^  and  a  Upschltz  continuous  Vg  with  ccmpact  support  such  that: 

Ug  -  Vg  or  B(zg,R);  (S(t)Ug)(z)  -  (s(t)vg)(z)  for  l*“*gl  <  1,  0  <  t  <  T|  and 

.  »  V?  .  for  KjAx,  kAy)  -z.  |<1,  0<n<N.  Applying  the  result  already  proved  (as 

j»k  3,K  0 

we  may  since  v",  S(t)Vg  also  have  compact  support.  Theorem  1  is  proved. 

Me  pause  to  comment  on  a  few  of  the  possible  extensions  of  the  preceding  results. 

First  of  all,  it  is  stral^tforward  to  treat  more  general  Hamiltonians  H(x,t,u,Du).  For 
exas^la,  in  one  space  dimension,  let  H(x,t,r)  be  I>ipshitz  continuous  in 
S  X  [0,T]  X  [-R,R]  for  each  T,R  >  0.  Then  an  approximation 


..n+1 


A.u" 


^ . 


G(x.,t  ,0 


J'  n'-J-p' 


is  consistent  if  g(x,t,a,. . . ,a)  ••  H(x,t,a)  and  monotone  if  is  a  nondecreaslng 

function  of  u"  , •  If  the  numerical  Masdltonlan  is  also  Lipschitz  continuous  on 
J~P  J’v' 

bounded  sets  of  R  x  [0,T]  x  we  can  again  estimate  d"  -  u(xyt^)  by  a  multiple 

1/2 

of  (At)  if  Ug  is  Lipschitz  continuous.  (The  simple  Propositions  1.1  and  3.1  need  to 
be  appropriately  generalized.  See  [11].) 

Next,  we  could  discuss  the  corresponding  stationary  problem  u  H(Du)  f(x)  in 
as  well  as  boundary  value  problew  (see  [9]),  but  we  will  not  formulate  any  precise 
results  hare.  It  is  also  clesr  that  implicit  approximations  can  be  handled  equally  well. 

Me  conclude  this  section  with  some  final  reswrks  in  the  context  of  the  equation 
u^  H(0a)  w  0  (trhich  apply  to  its  generalisations  as  well).  If  one  reexamin«s  the  above 
proofs  under  tlie  assusptlon  that  H  and  g  are  globally  Lioschitz  continuous,  one  sees 
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that  the  estiaate  on  -  u(Xj,yj^,nAt)  depends  on  Ufj  through  Ita  Lipechltz  constant 

L  *  L(Uq)  (provided  that  Uq  is  kept  bounded)  la  the  form 

l“j,k  ■  <  C{l.^(dt)’''^) 

(where  %ie  assuae  L  la  not  aaall  and  (At)  Is  not  large).  Using  this  fact  and  the 
nonexpanslve  nature  of  S(t)  (Proposition  1.1  (ID)  and  G  (Proposition  3.1  (111))  we 
conclude  tliat  If  U?  v"  ,  v  are  the  discrete  and  exact  solutions  for  Initial  data 

J#X  JfX 

Uq  (possibly  not  Lipachlts)  and  Vg  (Lipachlts),  then 

'“l.k  -  *  '"!.k  -  "'S.k'  ^  "'J.k  - 


♦  Iv(x^.y^.t^)  -  u(x^.y^,t„)l 
<  aiu.-v.l  ♦  C(L(w^))^(At)’''^  . 

U  0  O 


This  allows  us  to  conclude  the  convergence  of  the  nuaerleal  scheae  for  general 

Ug  e  BOCdf')  and  tha  approprlata  rata#  For  axMipXa«  If  la  H51dar  contlnuoua  with 

exponent  a,  we  can  choose  Vg  above  so  that  the  error  Is  at  aost  cons.  ( At)  ^*^*^"*** 
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Section  5.  Convergence  of  the  vanishing  viscosity  method. 

It  has  long  been  standard  to  attempt  to  approximate  (IV7)  by  the  problem 


(IVP), 


+  H(Du^)  -  edu*" 


0  in  R**  X  (0,“) 


u^(x,0)  -  Ug(x)  In  ***  . 

By  analogy  with  fluid  mechanics,  this  method  Is  referred  to  as  the  method  of  'vanishing 
1  ^  N  1  ^  M 

viscosity".  If  He  “o  ®  **  *  standard  results  and  methods  for  quasl- 

llnear  partial  differential  equations  yield  the  existence  and  uniqueness  of  a  solution  u^ 
of  (IVP)^  in  the  class  BUC(R^  *  [0,T]  )  f\  x  (0,T))  (l.e.,  continuous  second 

order  spatial  emd  first  order  time  derivatives)  for  all  T  <  >.  Our  main  result  Is 


Theorem  5.1;  Assume  H  is  locally  Llpschltz  continuous  on  Uq  Is  bounded  and 

Llpschltz  continuous  on  and  T  >  0.  Then,  if  u^  denotes  the  solution  of  (IVP)^ 

and  u  denotes  the  viscosity  solution  of  (IVP),  we  havet 


(5.1) 


sup  sup  |u®(x,t)  -  u(x,t)|  <  c  /e 


0<t<T  ^„N 
xSR 


where  c  depends  only  on  the  Llpschltz  constants  of  u^,  H  and  T. 


Proof;  To  simplify  the  presentation  we  will  only  consider  the  case  when  H(0)  ••  0  and  u, 
u^  -►  0  as  |x|  ♦  •  uniformly  in  t  6  (0,T)(T  <  »).  The  general  case  Is  easily  obtained 
by  (now)  routine  adaptations  of  this  simpler  case  ((1],  [2]). 

Assume 


(5.2) 


sup  sup  (u(x,t)  -  u  (x,t))  >  0 
xeR 


and  let  i|(  :  R**  x  r*  x  r**  x  r  be  given  liy 

(5.3)  i|((x,t,y,8)  “  u(x,t)  -  u®(y,s)  -  ^  (t+e)  +  (5M  ♦  |)8Jx-y,t-s) 

where  H  -  «  >  0  and  6^,  8  are  as  in  (4.7),  (4.9),  (4.21).  (We  now  use  x,  y  to 

denote  points  In  1^.)  Just  as  In  the  previous  case,  we  conclude  that  there  Is  a  point 
(Xg,tg,yg,8g)  which  maximizes  ♦  over  (R**  x  (o,T)  )^.  Moreover 
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(5.4) 


so  that 

(5.5)  •>‘o'yQ>^  *  ‘V®o'^  *  • 

As  before  we  will  consider  the  cases  tg,  Sg  >  0i  tg  ~  0,  Sg  >  0>  and  tg  >  0,  Sg  •  0 
separately.  Before  doing  so,  %ra  review  a  few  properties  of  u^.  We  have  the  elenentary 
estimates 

(5.6)  |u®|  <  M,  iDu'l  <  IDu„l  .  in  r"  *  (0.-) 

L  (r") 


(see,  e.g.,  [21,  [10] ).  Therefore  H(Du  )  is  bounded  Independently  of  e.  The  following 
lemma  will  then  allow  us  to  estimate  the  modulus  of  continuity  in  tijte  of  u^  in  the  form 
(5,7)  |u®(K,t)  -  u®(x,s)|  <  K  /e  lt-sl’‘^2+  Klt-s|  . 

Lemma  5.2.  Let  v  e  C^'’(l^*  x  (0,«))  A  m’'“(R*‘  x  (0,»))  satisfy  |v^  -  edvl  <  Kg  In 


R  X  (0,o).  Then  there  is  a  constant  K  depending  only  on  Kg  and  sup  IDvl 

t>0 


such 


L  (B  ) 


that  lv(x,t)  -  v(x,s)|  4  K(/e  |t-sl''^2+  |t-s|)  for  x  e  R** ,  t,s  >  0. 

^  11 

Proof.  Let  p  e  C  (R  )  be  a  standard  molllfler  supported  in  the  unit  ball  and  satisfying 

/  p(x)dx  -  1.  Put  v^  •  P„*v.  Clearly  lv_  -  edv^l  <  K-  in  R^.  Thus 
*  ad  dt  d  u 

'V'  -  H  ^  *^0  ^  ®  '*^a'  -  N 
L  (r")  L  (R  > 


+  ££ 

lov  1 

0 

d 

®  l’cr**) 

<  Ka 

+  S£ 

IDVI  _  „ 

0 

d 

•  M 

I.  (e") 

for  t  >  0 


where  c  depends  only  on  p.  Therefore 


fv(x,t)  -  v(x,s)|  <  |v(x,t)  -  V  (x,t)|  +  |v(x,s)  -  v  (x,8)|  +  (K-  +  ^)|t-s| 

d  d  0  d 


<  K(o  +  |t-s|  ♦  ■!  |t-s|)  , 

where  K  denotes  several  constants  with  allowed  dependencies.  Setting  a  <•  |t-n|^ 

yields  t)ie  result. 

We  now  turn  to  the  eases  tg  >  0,  Sg  ■  0i  and  tg  •  0,  sg  >  0.  Here  only  continuity 
considerations  are  Involved,  as  before.  If  tg  >  0  and  Sg  ■  0  we  have 
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<• 


5M  +  O  <  i()(Xg,tg,y|j,0)  <  5M  +  7  +  |u(y„,t^)  -  u(K„,t„) 


0'  0' 


0'  0 


+  |u(Xq,0)  -  u{Xg,tu)| 


<  5M  +  -  ♦  c(|Xj|-yg|  +  t|j)  . 


Now,  In  a  manner  similar  to  (but  simpler  than)  the  3rd  case  In  Section  4,  we  can  conclude 

2  2 

that  |xg~ygl  4nd  tg  are  bounded  by  multiples  of  a  and  so  o  <  ca  In  this  case.  If 
Vj, 

a  ~  e  * ,  we  have  the  desired  estimate.  The  case  tg  •  Sg  -  0  Is  subsumed  under  the  one 
just  treated.  If  tg  -  0  and  Sg  >  0,  we  have 

"  ‘5M+<l/2)e^(XQ-yj,,-Sg)  +  ^  •g  < 

u®(yg,0)  -  (5M+o/2)B^(X(j-yg,0)  . 

From  this  and  (5.5),  (5.7)  we  deduce  an  estimate 

.2 


-2-  <  K(/^  /T  +  8„) 


where  K  Is  Independent  of  a.  This  Implies  tl»at  Sg  <  ♦  a^).  Now  5M+o  < 

♦  (Xg,0,yg,Sg)  4  5M  ♦  ^  ♦  Ko^  +  K/e  /s^  +  KSg.  Using  the  previous  est^Jiate  and  letting 


e/4 


4  Ke''2  . 


a  ••  e'"  we  conclude  that  a 

The  final  case,  Sg  >  0,  tg  >  0,  usee,  the  equations.  From  the  fact  that  u  is 
Llpschltx  continuous  In  x  and  t  we  deduce  that  ^  ^ ' 

|d^0  (x.-y.,t„-S-)  I  are  bounded  where  D_,  D,.  refer  to  the  derivatives  of  B  (x,t).  This 

t^ttUUUO  XV  Cl 

Implies 


'V^o''  'V*o'  * 


(5.8) 

Using  that  u  Is  a  viscosity  solution  we  have 


(5.9) 


--  (5MtC/2)(D^B^)(Xg-yg,tg-Sg)  ^ 


H(-(5M+C/2)(D^6^)(Xg-yg,tg-8g))  4  0 


€  2  1 

On  the  other  hand  using  u  8  C  '  and  that  (^g'^g)  minimises  (y,s)  *  4i(Xg ,tg ,y,s) 
over  X  [o,T]  we  also  have 
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<5.10) 


% 


<• 


(5.11)  Du®(yj,»j|)  +  {5M+o/2)(D^B|j)(Xjj-y|j,t^-8|j)  -  0  , 

and 

(5.12)  Au®(yg.ag)  -  (5M1-0/2) ( 4$^) (Xg-yj,tQ-ag )  >  0  . 

c  c  c 

Now  procMd  by  using  u^  +  H(Du  )  -  cAu  >  0,  and  (S.10),  (5.12)  to  daduce  that 
e(5(tto/2)<AB^)(Xj,-yQ,tj,-Sg))  -  H(-(5M+o/2)(D^B^)(X|j-yu,tp-8jj)) 

+  15+  (5M><l/2)(D^^B^)(Xg-y5,tQ-»o))  <  0  . 

Now,  using  (5.9),  ws  conclude  that 

and  this  yields  a  <  KC/a^.  Again  it  a  -  we  have  the  desired  estimate. 

RemarK.  By  contrast  with  the  analysis  in  section  4,  the  Upschits  continuity  of  H  was 
used  only  to  assert  that  (IVP)^  has  a  saxioth  solution,  if  H  is  Mrely  continuous  one 
still  has  tt*  and  Au*  in  x  (O,*))  for  1  <  p  <  •  and  the  estimate  (5.1)  can 

still  be  proved. 
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M.  ABSTRACT  fCaallmia  an  ravaraa  dda  II  naaaaaaiF  and  Idinlllir  hr  Uaek  mamhar) 

Equations  of  Hamilton-Jacobi  type  arise  in  many  areas  of  application, 
including  the  calculus  of  variations,  control  theory  amd  differential  games.  The 
associated  initial-value  problems  almost  never  have  global-time  classical 
solutions,  and  one  must  deal  with  suitable  generalized  solutions.  The  correct 
class  of  generalized  solutions  has  only  recently  been  established  by  the  authors. 
This  article  establishes  the  convergence  of  a  class  of  difference  approximations 
to  these  solutions  by  obtaining  explicit  error  estimates.  Analogous  results 
are  proved  by  similar  means  for  the  method  of  vanishing  viscosity.  I 
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